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ABSTRACT 

Context. Chemically peculiar stars of the upper part of the main sequence show periodical variability in line intensities and continua, 
modulated by the stellar rotation, which is attributed to the existence of chemical spots on the surface of these stars. The flux variability 
is caused by the changing redistribution rate of the radiative flux predominantly from the short-wavelength part of the spectra to the 
long-wavelength part, which is a result of abundance anomalies. Many details of this process are still unknown. 

Aims. We study the nature of the multi-spectral variability of one of the brightest chemically peculiar stars, 0 Aur. 

Methods. We predict the flux variability of 6 Aur from the emerging intensities calculated for individual surface elements of the 
star taking into account horizontal variation of chemical composition. The surface chemical composition was derived from Doppler 
abundance maps. 

Results. The simulated optical variability in the Strbmgren photometric system and the ultraviolet flux variability agree well with 
observations. The lUE flux distribution is reproduced in great detail by our models in the near ultraviolet region. A minor disagreement 
remains between the observed and predicted fluxes in the far ultraviolet region. The resonance lines of magnesium and possibly also 
some lines of silicon are relatively weak in the ultraviolet domain, which indicates non-negligible vertical abundance gradients in the 
atmosphere. We also derive a new period of the star, P = 3.618 664(10) d, from all available photometric and magnetic measurements 
and show that the observed rotational period is constant over decades. 

Conclusions. The ultraviolet and visual variability of 6 Aur is mostly caused by silicon bound-free absorption and chromium and iron 
line absorption. Manganese also contributes to the variability, but to a lesser extent. These elements redistribute the flux mainly from 
the far-ultraviolet region to the near-ultraviolet and optical regions in the surface abundance spots. The light variability is modulated 
by the stellar rotation. The ultraviolet domain is key for understanding the properties of chemically peculiar stars. It provides a detailed 
test for surface abundance models and comprises many lines that originate from states with a low excitation potential, which enable 
detecting vertical abundance gradients. 
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1. Introduction 

The atmospheres of main-sequence A and B stars are strongly 
affected by the radiative force. The radiative force acts selec¬ 
tively on some heavier e lements, which leads to atomic diffusion 
and atomic segregation (IVauclairil2003t |Michaudll2004l) . In nor¬ 
mal AB stars the atomic diffusion is counteracte d by the winds 
in ho tter stars and by turbulence in cooler stars dMichaud et alJ 
|201 111 , consequently, their abundance anomalies are only mod¬ 
est. The atmospheres of chemically peculiar (CP) stars are 
more stable, however, and consequently show distinct chemi- 
cal composition remarkably different from the normal one (see 
lRomanvukl2Q07[ for a review). 

The magnetic field probably is the stabilizing force in a 
group of CP stars called magnetic chemically peculiar stars. In 
these stars the magnetic field of the order of hundreds of Gauss 
or higher dominates the whole atmosphere. Normal A stars may 
possibly also have a surface magnetic field, but with an intensity 
lower by two orders of magnitudes, as s hown for Vega and Sirius 
dLignieres et alJ2dO^IPetit et al.ll2Qll1) . Such weak fields do not 
significantly affect the atmospheres. There seems to be a clear 
dichotomy between stars with strong and weak magnetic fields 
dAuriere et alJ120071) . Consequently, it seems that the so-called 
normal AB stars and the peculiar ones are in many respects very 
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similar apart from the strength of their anomalies (the strength of 
the magnetic field and the difference in chemical composition). 


This view may be supported by the observation of the light 
variability of the A star. The light variability of CP stars is caused 


by the flux redistribution due to bound-free (ionization. Peterson 

1970 

; iLanz et al. 19961) and boud-bound (line. IWolff & Wolfi 

1971 

; iTrasco 19721 Molnar||1973|) transitions in surface abun- 


dance spots, which is modulated with stellar rotation. It is tempt¬ 
ing to attrib ute the milli-m agnitude variability observed in nor¬ 
mal A stars dBalonal2013l) to the same mechanism as in CP stars. 


However, the test of this hypothesis is probably beyond the 
possibilities of current observational techniques. The situation 
in CP stars is much simpler. The theoretical light curves can be 
predicted using model atmospheres calculated with actual abun¬ 
dances from abundance maps of individual stars. A comparison 
of such theoretical light curves with observed one s supports the 
current picture of the light variability of CP sta rs dKrticka et al.l 
120091: IShulvak et al.l l2()l^ iKrticka et al.ll20l3) . This approach 
cannot be used for rotationally variable normal A stars be¬ 
cause the expected abundance and line profile variations are pro¬ 
hibitively small for abundance mapping. 

Consequently, studies of the light variability of AB stars 
should concentrate on CP stars with amplitudes of the order of 
hundredths of magnitudes. But our knowledge is still far from 
complete even in this field. Most studies concentrate on hotter 
CP stars, s UMa being the only example of a cooler CP star stud- 
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Table 1. List of the lUE observations of 0 Aur 


Image 

SWP camera 
Julian date 
2,400,000 -h 

Phase 

Image 

LWR camera 
Julian date 
2,400,000 -h 

Phase 

21327 

45630.77975 

0.068 

2111 

45630.78658 

0.070 

21331 

45631.04433 

0.141 

2135 

45632.79204 

0.624 

21332 

45631.07042 

0.148 

2139 

45633.03122 

0.690 

21353 

45632.79787 

0.625 

2150 

45634.82937 

0.187 

21357 

45633.03467 

0.691 

2154 

45635.06043 

0.251 

21358 

45633.06212 

0.698 

2166 

45636.79184 

0.729 

21375 

45634.79290 

0.177 

2184 

45638.78399 

0.280 

21378 

45635.06294 

0.251 

2186 

45638.89943 

0.312 

21390 

45636.78577 

0.727 




21412 

45638.77651 

0.278 




21414 

45638.90422 

0.313 





Notes. Phases were calculated according to Eq. (|6j. 


ied so far (IShulvak et al.ll201^ based on maps of lLiiftinger et alJ 
l2003h . Moreover, because we lack observational data in the UV, 
it is not always possible to couple the study of optical light vari¬ 
ability with that of the UV spectral energy distribution variabil¬ 
ity. However, this is crucial for CP stars, in which the light vari¬ 
ability originates in the UV, and the optical variations are just a 
glimpse of it. 

To fill this gap, we here provide a study of the optical 
and UV flux variability of one of the brightest CP stars, 0Aur 
(HR 2095, HD 40312). This star belongs to a group of well- 
studied CP stars. The Balmer line spectroscopy revealed the in- 
fluence of the Lorent z force on the stellar atmosphere of this star 
dShulvak et alJ l2007h. The star is a fre q uent target of Doppler 
mapping (e.g., fpiokhlova et alJ [T986I: iRice & Wehlaul 1 199 it 
lKuschniall998b ’ lRice et al* 2004i) . We make use of the surface 
abunda nce maps of variou s elements derived using Doppler im¬ 
ages bv iKuschnid (Il998ah and simulate the flux variability using 
model atmospheres. 


2. Observations 


g Aur has been recognized as an AOp Si star by ICowle y et al. 
(19^). The variability in UBV was established by Iw inzen 
(11974 ) and confirm ed by two sets of uvby measu r ement s 
of AdelmanI (Il997h and Adelman & KaewkornmaunS (l2005h . 
The star w as also observed by Hipparcos in the Hp band 
(lESAlll9^). The measurements of t h e effective magnetic held 
made by Landst reet & Bo^ (Il977h : IWade et aP (120001) . and 
[Silvester et alJ (l2012l) renre^nt further valuable phase informa¬ 
tion on the variability of 6 Aur. 

We supplemente d these data by lUE sp ectra extracted from 
the INES database (I Wamsteker et al.l l2000l see Table [1]) using 
the SPLAT package. We selected mostly high-dispersion large- 


aperture spectra in the domains 1150-1900 A (SWP camera) and 
2000-3000 A (LWR camera). We used phase-dependent fluxes 
in our analysis. Moreover, we calculated artificial magnitudes 
from the spectra by convolving them with a Gaussian Alter with a 
dispersion of 25 A (for SWP spectra) and 100 A (for LWR spec¬ 
tra). 


3. Phenomenological modelling of the light and 
magnetic variability 

The aim of the phenomenological modelling of the 0 Aur light 
and magnetic held variability is to improve the rotational period 



Fig. 1. Light curve of 0 Aur in Hp fitted by the model light curve 
(full line), which is the sum of two partial light curves corre¬ 
sponding to individual spots (dashed lines) centred at tpoi and 
ipo 2 - Observations are denoted by circles. 


of the star using all relevant sources of phase information, the 
few parametric descriptions of the light-curve shapes in the spec¬ 
tral domain from 1300 to 5500 A, and to determine the magnetic 
held geometry. 

3.1. Phenomenological model 

We found that each of the studied 19 light curves with effec¬ 
tive wavelengths from 1300 to 5500 A is a smooth single wave 
curve, which can be approximated by a harmonic function of 
second order. Nevertheless, the shapes of the light curves in dif¬ 
ferent bands are apparently different. We applied the weighted 
advanced principle component analysis (APCA) to the param¬ 
eters of the harmonic fits, which allowed us to And hidden re¬ 
lationships among th e individual light curves (for details see a 
brief introduction in lMikulasekll2007h . We concluded that all 
light curves can well be expressed by a linear combination of 
only two principle components and that the amplitudes of the 
third and higher principal components are negligible. 

This allowed us to build a relatively simple two-component 
phenomenological model with a minimum of free parameters 
valid for all studied light curves, where each of them is expressed 
by a linear combination of two principle normalized symmetrical 
functions (see Eig.[T]i. The models of the light curves are 


F(d,A)=Aok-2.5 logioU, 


where 


7i - 



1 - cosh 



2.288 d 


( 1 ) 

( 2 ) 


The phase models are defined as 

^ = ^ ; (fj = (d-(fioj)-round(d-(fioj), ; = 1, 2, (3) 


where t? is the so-called phase functio n (sum of the commo n 
phase ip and the epoch £, for details see iMikulasek et aD2008l) . 
and f is the time in BJD. The origin of the phase function (zero 
phase time) Mq was fix ed at 2450001.88 1 (the basic magnetic 
minimum according to I Wade et al.ll2000h . P is the period (as¬ 
sumed to be constant), ipoi and ^02 are phases of centres of basic 
phase profiles (the phases of photocentres of photometric spots), 
d\ and d 2 are the parameters expressing the width of the model 
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Table 2. List of the phase data sources used in the phenomeno¬ 
logical modelling of 0 Aur variability. is the mean epoch ac¬ 
cording to our ephemeris (see Eq.j^. N is the number of obser¬ 
vations. 


Year 


N 

Em 

O-C 

Source 

1972 

UBV 

45 

-Till 

-0.07(7) 

IWin 7 eiUl 974 ) 

1977 

Si 

19 

-1931 

0.07(6) 

iLandstreet & Borral n977h 

1984 

lUE 

106 

-1207 

-0.05(10) 

this paper 

1992 

Hp 

131 

-398 

0.03(3) 

Ies^IT^ 

1994 

uvby 

188 

-210 

- 0 . 02 ( 2 ) 

lAdelmanll1997f 

1998 

C 

11 

209 

-0.07(4) 

IWade et all 12000h 

2004 

uvby 

276 

859 

- 0 . 01 ( 2 ) 

Adelman et al. (2005) 

2008 

Si 

7 

1206 

0 . 00 ( 2 ) 

Silvester et al. (2012) 


spot profiles in phases, Ajx are the amplitudes of the light vari¬ 
ations for the y'-th spot, (fj is the phase counted from the centre 
of the j-th profile, and Aoi is the mean magnitude of the k-i\\ 
set of observations from an individual author and an individual 
photometric colour. 

The simplest phenomenological model of the phase curve 
of the longitudinal magnetic field (corresponding to the oblique 
central dipole) is the cosinusoid with the extrema at the phases 
when the magnetic poles are passing the stellar meridian. 
Nevertheless, the recent observations of the longitudinal mag¬ 
net! cfield_B|_of0_Aurindicate_amorecomplex phase behaviour 
(see lWade et alJ20()(itlSilvester et al.l2012L and also Fig. |5]l. That 
is why we approximated the magnetic phase curve by the har¬ 
monic polynomial of the second order in the following special 
form, which allows us to estimate the extent of the distortion of 
the basic purely cosine course: 

L(i?) - B\+ Aj, cos (2 TT Bi) + A 4 cos (4 tt d-z) + (4) 

-I- A 5 [sin(2 n- 5 sin(4 n )? 3 )j, §3 - - (po 3 . 


where ^03 is the phase of the magnetic curve minimum, Bi is 
the mean longitudinal magnetic field, A 3 is the semi-amplitude 
of the cosine component of the magnetic variations. The semi¬ 
amplitudes A 4 andAs quantify the significance of the possible 
symmetrical/antisymmetrical distortions of the observed phase 
curve. 

The model describes all observed light curves obtained in 
18 different photometric colours and magnetic variations with 
a fairly high fidelity with only 67 free parameters. This phe¬ 
nomenological model was applied to observational data repre¬ 
senting 771 individual measurements. The measurements more 
or less evenly cover the time interval 1971 - 2008 (see Table |2]), 
and consist of three sets of magnetic measurements, UBV ob¬ 
servations, two sets of uvby measurements, Hp photometry, and 
106 magnitudes in 11 ultraviolet wavelengths derived from lUE 
spectra (see Sect. 0 . 

We determined individual parameters a and their uncertain¬ 
ties of the model described by Eqs. ([T]i and (|4li by deriving for 
each measurement its model prediction F{ti, a) using the stan¬ 
dard weighted least-square method applied to the set of mea¬ 
sured values (magnitude, magnetic field) y,, with the estimated 
uncertainty cr,-. 


Ay,- = y,- - F{ti, a); 



Va" = 0; 



For details see also lMikulasek et ^ (1201 lah . 


(5) 



Years 


Fig. 2 . O-C diagram of 0Aur indicates that its period is 
nearly constant. Circles indicate 0-Cs derived from photome¬ 
try, squares correspond to 0-Cs found from magnetic measure¬ 
ments. 


3.2. New ephemeris and O-C diagram 

One of the principle results of our phenomenological modelling 
was the improved period of the observed stellar variations: P = 
3‘?618 664(10). The time of the zero phase and the zero epoch, 
Mq - 2 450001.88 1 , were adopted from the ephemeris pub¬ 
lished in IWade et al.l (l2000ll . The times of the zero phase 0{E) 
can be approximated by the relation 

0(E) = 2 450 001.881 H- 3^618 664( 10) x £, ( 6 ) 


where E is the integer, also called the epoch. The pe¬ 
riod agrees well with other determ inations (IWade et alJ I 2 OOOI: 
lAdelman & Kaewkornmaun3l2005h. The phase sh ift between 
this ephemeris and that used by iKuschnia (Il998ah in the time 


of their spectral observations is Atp — -0.268. 

We also tested a possible period variability. We assumed that 
the phase function is approximated by the relation 




t-Mp 

Po 


P 

2 


-Mq Y 

Po I 


0 - Mq + PqE + 


PPoE^ 


(7) 


where Pq is the instant period at Mq and P is the time derivative 
of the period. We found that P - 3(5) x 10“®, so we can conclude 
that there is zero quadratic term in the 6 Aur ephemeris and the 
period is constant. 

A more complex behaviour of the period can be revealed by 
means of the O-C diagram inspection. The phase function can 
be approximated by the relation 


?7(f, r) 


t — Mo — (0-C)r 


( 8 ) 


where Mq and Pq are fixed, Pq - 3‘?618 664, and (0-(J)r are the 
mean O-C v alues of the subsets r of the observational data (for 
details see in iMikulasek et al.l201 Ibh . The O-C diagram in Fig.|^ 
(constructed using observations in Table|2]i does not indicate any 
period variations during the past 33 years. 


3.3. Light curves 

We present the light curves of 0 Aur in the region 1300 - 5500 A, 
which covers a substantial part of the stellar flux, in Fig.|3] The 
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Fig. 3. Observed light curves of 0 Aur in UV {upper panel) and 
in standard photometry {lower panel). Symbols on the right side 
denote the effective wavelengths in A or a hlter. Point areas in 
the lower panel correspond to the weight of measurements. Grey 
lines are the results of phenomenological modelling. 



Fig. 4. Amplitudes of two photometric spots centred at 1^02 = 
0.11 (squares) and ipoi - 0.39 (full dots) at various wavelengths. 


Fig. 5. Variations of the longitudinal component of t he magnetic 
held o f 0Aur. Squares denote the observation of Wade_eL^ 
( I 2 OOOI) . diamonds the observati on of Landstreet & Borral(ll977h . 
and circles the observation of [Silvester et alJ (l2012l) . The grey 
line is the symmetric model ht. 


largest amplitudes are shown by the light curves at the short- 
wavelength end of the spectral region. The phenomenological 
modelling of the light curve set by the linear combination of 
two basic symmetric phase prohles (see Eq.[T]| was successful, 
the model is able to plausibly ht the shapes of all observed light 
curves. The dominant feature of the spectrum is centred at the 
phase ^01 = 0.389(4). Our test showed that the width d 2 of the 
second prohle centred at (fo 2 = 0.112(10) is nearly identical to 
the width of the hrst prohle. Consequently, we assumed that the 
widths of all phase structures are the same, d\ - d2 - d - 
0.165(3). 

The observed behaviour of light curves can be interpreted by 
the presence of two dominant photometric spots on the surface 
of the rotating star whose photo-centres pass the stellar median 
at phases 0.389 and 0.110. The contrasts of the spots with respect 
to their surroundings and the contrast of the spots themselves de¬ 
termine the shape of the light curve in a particular photometric 
colour. The principle spot is bright only for wavelengths longer 
than 3400 A, otherwise it is dark. The second spot is bright for 
wavelengths longer than 1850 A. For shorter wavelengths, both 
spots are dark, and together they create one extensive shady fea¬ 
ture. 


3.4. Geometry of the magnetic field and its relation to 
photometric spots 

The magnetic held variations are described by the model Eq. dUi. 
The phase of the sharp magnetic minimum is <^03 = 0.024(4). 
The mean value of the longitudinal component of the magnetic 
induction of 0 Aur is positive (Z?i = 85(4) G), which means that 
the positive pole of the magnetic dipole is in the northern hemi¬ 
sphere of the star. The derived semi-amplitude of magnetic vari¬ 
ations is A 3 = -243(7) G, and the semi-amplitude of the sym¬ 
metrical component A 4 = -44(5) G is nonzero, contrary to the 
semi-amplitude of the antisymmetric term A 5 = -8(10) G. This 
means that the phase curve of the magnetic variations is deh- 
nitely distorted, but symmetric around the phase of the magnetic 
minimum (see Fig.|5]l. This form of the phase curve can be ex- 
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plained by assuming that the magnetic dipole is slightly shifted 
from the centre in the direction of the negative magnetic pole. 

The phases of the centres of assumed photometric spots nei¬ 
ther coincide with the phases of magnetic held extrema nor with 
phases when the magnetic equator passes the centre of the vis¬ 
ible stellar disc. More over, the spots wit h the highest element 
abundances mapped by iKuschnig (Il998all appear on the visible 
disc around the phases 0.3 - 0.4, again showing a mismatch with 
the surface magnetic held. Thus we conclude that there is no 
obvious connection between the photometric structures on the 
stellar surface and the stellar magnetic held geometry. 


4. Calculating the flux variability 

4.1. Model atmospheres and synthetic spectra 

The stellar parameters of 6 Aur together with the ran ge of sur¬ 
face ab undances of individual elements taken from lKuschni3 
(Il998al) are given in Table [3] The effective temperature and sur¬ 
face gravity were derived from htting the Hb lines. The pro¬ 
jected rotational velocity was derived from spectroscopy. The 
range of surface abundances was extracted from surface abun¬ 
dance maps, which were derived by Doppler imaging. In our 
study the abundances are dehned relative to hydrogen, that is, 

fiel = log(Ael/AH). 

We used the code ATLAS 12 for the model atmosphere 
calculations (iKurucj l2005t ICastellil l2005l) . For each model 
atmosphere we calculated the synth etic spectrum us i ng th e 
code SYNSPEC with atomic data from iLanz & Hubenvl(l2007l) . 
The adopted SYNSPEC line list was enlarged ^ the data 
for chromium and iron taken from Kurucz (2009U We com¬ 
puted angle-dependent emerging intensities for 20 equidistantly 
spaced values of p - cos 6 , where 0 is the angle between the 
normal to the surface and the line of sight. 

In our models we assumed fixed effective temperature and 
surface gravity (see Table [3]l and adopted a generic value of the 
microturbulent velocity Vturb = 2kms“*. The abundance of Mg, 
Si, Cr, Mn, and Ee was selected according to the surface abun¬ 
dance maps (s ee below). We assum ed solar abundance of the 
other elements dAsnlund et al.ll206^ . 

The model atmospheres and the angle-dependent inten¬ 
sities /(/1,0, fiMg, esi, ecr, fiMn, epe) Were calculated for a five- 
parametric grid of magnesium, silicon, chromium, manganese, 
and iron abundances (given in Table [Hi. Eor silicon and iron 
this grid fully covers the range of abundances found on the sur¬ 
face of 6 Aur (Table [3]). Eor magnesium, chromium, and man¬ 
ganese the lowest abundances were omitted from the grid be¬ 
cause these elements with such low abundance do not influence 
the emerging fluxes. This helped us to signihcantly reduce the 
required computer time. Our test showed that the highest abun¬ 
dances of helium, titanium, and strontium on the surface of 0 Aur 
(see Table O are so low that they do not influence the emerging 
flux. Consequently, these elements were excluded from the grid 
of calculated abundances in TableH] and we only assumed a sur¬ 
face mean abundance of these elements. 


4.2. Phase-dependent flux distribution 

The derived grid of spectra (Table |4]i was used to calculate the 
radiative flux in a bandpass c at the distance D from the star 

' http://kurucz.harvard.edu 


Table 3. 0 Aur parameters derived bv lKuschnid(ll998ah . 


Effective temperature 

10 500 K 

Surface gravity log g (cgs) 

3.5 

Inclination i 

60° 

Rotational velocity projection Vrot sin i 

55 km s“* 

Stellar radius* R, 

4.34(5) Ro 

Helium abundance 

-2.8 < SHe < -1.9 

Magnesium abundance 

-7.5 < fiMg < -3.2 

Silicon abundance 

-3.6 < Ssi < -2.8 

Titanium abundance 

-7.9 < St, < -6.9 

Chromium abundance 

-6.0 < Scr < -3.4 

Manganese abundance 

-6.2 < SMn < -4.7 

Iron abundance 

—4.2 < £rpe < —3.2 

Strontium abundance 

-9.1 <esr<-7.7 


Notes. Derived in Sect.[7l 

Table 4. Individual abundances SMg, esi. scr, eMn. and epe of the 
model grid 


Mg 

-4.1 

-3.6 

-3.1 

Si 

-3.7 

-3.2 

-2.7 

Cr 

-5.9 

-5.4 

-4.9 

Mn 

-5.7 

-5.2 

-4.7 

Fe 

-4.2 

-3.7 

-3.2 


Table 5. Coefficients of analytical hts of the photonic passbands 
of the Stromgren photometric system Alters uvby Eqs. (fTTb and 



fli 

a2 

^3 

^4 

df, 

06 

dc 

CTc 

u 

21.3 

-14.3 

3590 

8.85 

448 

-1430 

3475 

10" 

V 

94.4 

-652 

2000 

55.2 

305 

0 

4100 

10" 

b 

76.7 

-268 

0 

111 

-896 

2790 

4675 

10" 

y 

24.8 

465 

0 

47.1 

578 

0 

5475 

10" 


(lMihalaslll978h 

f 4(e,Q)cos0dD. (9) 

visible 

surface 

The intensity Ic(6,0.) at angle 0 with respect to the normal 
at the surface was obtained at each surface point with spher¬ 
ical coordinates Q by interpolating between the intensities 
hfG, SMg, ssi. ecr, SMn, Spe) Calculated from the grid as 


Ic{Q, SMg, esi, ecr, eMn, epe) = 

<t)c(/l) I{A, 0, fiMg, esi, ecr, SMn, Spe) dd. (10) 

The response function <[)c(/l) of a given bandpass c is a Gauss 
function for the UV magnitudes (the same as for lUE fluxes, see 
Sect. |2]i. Eor the Stromgren photo metric system w e use analyti¬ 
cal fits to photonic passbands from iBesselll (1201 ll) . 



0 ,( 4 ) = 


exp{—x[aic + x(a 2 ca^icx)]}, x < dc, 
exp{—.r [a4c -P xia^c + a(,cxy\], x > dc. 


( 11 ) 


where 


X - 


(A - dcf 


erf. 


( 12 ) 


where A is the wavelength in A. The coefficients of the fits are 
given in Table |5] 
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"^ross 


Fig. 6. Upper plot: Temperature dependence on the Rosseland 
optical depth Tross in the reference model atmosphere with eMg = 
-4.1, ssi = -3.7, ficr = -4.4, fiMn = -5.7, and spe = -4.2. 
Lower plot. Difference between temperature in the model atmo¬ 
spheres with abundance of the selected element enhanced by a 
factor of 10 and the temperature in the reference model atmo¬ 
sphere. 



wavelength X [A] 

Fig. 8. Magnitude difference between the emerging fluxes 
calculated with an enhanced abundance (by a factor of 10) of a 
selected element and the reference flux (Ea.[T4li. The fluxes 
were smoothed by a Gaussian Alter with a dispersion of 100 A. 


be seen from the plot of the magnitude difference 


Ami = -2.5 log 


77i(^>Mg7 ^Fe) 


trf 


(14) 


The magnitude difference between the flux fc at a given 
phase and the reference flux in bandpass c is defined as 

Am, = -2.5 log (Aj. (13) 

The reference flux is obtained under the condition that the mean 
magnitude difference over the rotational period is zero. 

5. Influence of heavier elements on atmospheric 
structure 

Heavier elements significantly affect the structure of stellar at¬ 
mospheres. Although their abundance is much smaller than that 
of hydrogen even in CP stars, their contribution to opacity is 
crucial. Bound-free transition of silicon and numerous line tran¬ 
sitions of iron group elements (in our case, chromium and iron) 
are able to absorb the stellar radiation, which leads to heating 
of continuum-forming regions (Tross ~ 0.1 - 1). This is shown 
in Fig. 12 where we plot the temperature difference between the 
models with enhanced abundance of a selected element by a fac¬ 
tor of 10 and a reference one. 

The heavier elements also lead to a redistribution of the 
emerging flux typically from the far UV to the near UV re¬ 
gion and to the optical domain. This is shown in Fig. |7] where 
we compare the emerging fluxes calculated in a model with en¬ 
hanced abundance of the selected element by a factor of 10 and 
the reference one. Silicon redistributes the far UV flux to the 
spectral regions with A > 1700 A. The influence of chromium 
and iron on the emerging flux is more complex, but in general, 
these elements redistribute the flux from the far UV region to the 
spectral region with A > 2900 A. Magnesium and manganese do 
not affect the emergent flux significantly. 

The flux changes caused by the heavier elements manifest 
themselves in the changes of the apparent magnitude. This can 


in Fig. m Hf is the reference flux calculated for eMg = -4.1, 
esi = -3.7, ecr = -4.4, cmh = -5.7, and ep, = -4.2. For sil¬ 
icon the magnitude difference decreases with increasing wave¬ 
length. There is a minimum of Am^ at about 5200 A for iron 
and at 5400-5600 A for chromium, supporting the conclusion 
that these elements signifl cantly contribute to the flux depres¬ 
sion at these wavelength s (iKupka et alj|2003l iKhan & Shulvald 
120071 : l&ticka et al.]l2009l) . As follows already from Fig. |7l mag¬ 
nesium and manganese do not signiflcantly affect the emerging 
fluxes in the case studied here. 


6. Predicted broad-band variations 

Different parts of the stellar surface are exposed to the observer 
as a result of the stellar rotation, which causes the stellar light 
variability. To account for the stellar rotation, predicted light 
curves were calculated from Fq. [13] as a function of the rota¬ 
tional phase. 

In Fig. 12 we provide the light variations calculated from the 
abundance map of one element alone assuming a fixed abun¬ 
dance of other elements (eMg = -4.1, esi = -3.7, ecr = -5.9, 
SMn = -5.7, and ep, = -4.2). From Fig.j^it follows that mainly 
silicon, chromium, iron, and partly manganese cause the light 
variability. These elements are able to affect the emerging flux, 
and there is a large overabundance of these elements in the spots. 
The light curves of individual elements reflect the variations of 
line equivalent widths. As a result of the different surface distri¬ 
bution of individual elements, the line equivalent widths curves 
and also individual light curves in Fig.|2differ in shape. 

Silicon bound-free absorption dominates the flux distribu¬ 
tion below roughly 1700 A. The absorbed flux is redistributed to 
longer wavelengths, consequently, the silicon flux variations at 
1300 A and 2500 A are in anti-phase in Fig. |2 For similar rea¬ 
sons, the UV flux variations and the visual light variations are in 
anthi-phase for chromium, manganese, and iron. 
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Fig. 7. Upper plot: Emerging flux from a reference model atmosphere (sMg = -4.1, esi = -3.7, sq = -4.4, SMn = -5.7, and 
spe = -4.2). Lower plot: Difference between the emerging flux from the model atmospheres with higher abundance of the selected 
element by a factor of 10 and the flux from a reference model. All fluxes were smoothed by a Gaussian filter with a dispersion of 
10 A to show the changes in continuum. The passbands of the uvby photometric system are also denoted. 




Fig. 11. Comparison of the predicted (solid line) and observed (dots) UV light variations for different wavelengths (marked in the 
graph). The light curves were derived from lUE spectra using a Gaussian filter with a dispersion of 25 A (for SWP spectra) and 
100 A (LWR spectra). Curves for individual wavelengths were vertically shifted to better show the variability. 


Taking into account the surface distribution of all significant 
elements (see Table |4]i, the observed and predicted light curves 
in the b and y filters of the Stromgren photometric system agree 
well in Eig. [10] The observed curves have slightly lower am¬ 
plitudes than the predicted curves in the v filter, whereas there 
is also a slight phase shift between these curves in the u fil¬ 
ter. A similar behaviour also appears in the UV light curves in 
Eig. [TT] Here we plot the observed and predicted light curves, 
which were processed using a Gaussian filter. The observed and 


predicted light curves agree very well in the long-wavelength 
part of the UV spectrum with A > 1600 A, whereas the observed 
curves have lower amplitudes than the predicted curves in the 
far-UV domain with A < 1500 A. 

The inhomogeneous surface distribution of individual ele¬ 
ments causes bright spots on the stellar surface. The spots, whose 
surface distribution can be derived using abundance maps and 
model atmospheres (see Eig. fTZb . cause the light variability. 
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Fig. 9. Predicted UV and visual flux variations of 0Aur calcu¬ 
lated using abundance maps of one element alone. The abun¬ 
dance of other elements was fixed. Light curves for selected 
wavelenghts in the UV domain (denoted in the graph) were cal¬ 
culated using a Gaussian filter with a dispersion of 100 A. The 
light variations were calculated for the u and y filters of the 
Stromgren photometric system. Light curves in individual bands 
were vertically shifted to better show the variability. 


Given the relative proximity of the star, 0Aur can be an ide al 
candidate for interferometric analysis (cf. lShulvak et al.l[2014l) . 



Fig. 10. Predicted light variations of 0Aur (solid lines) in the 
Stromgren photometric system computed taking the surface dis¬ 
tribution of all elements. The observed light variat ions (dots) 
are taken from I Adelman & Kaewkornmaungl (l2005h . The light 
curves in individual filters were vertically shifted to better show 
the light variability. 




Fig. 12. Emerging intensity from individual surface elements of 
0 Aur at various rotational phases. Upper panel; visible y band. 
Lower panel: UV band centred at 2700 A. Both for fj. - 1. 


7. Ultraviolet flux variations 

We have shown that most of the light variability of 0Aur can 
be explained by the flux redistribution from far UV to the near 
UV and visible regions. This is supported by the comparison of 
UV and visual light curves, which are typically in anti-phase. 
However, there is some disagreement between observed and 
predicted light curves. The predicted M-amplitude is too low 
(Fig . [Tol l and, in addition, the predicted variations in the far-UV 
are stronger than the observed variations (Fig. [IB- Further dis¬ 
cussion of ultraviolet fluxes is necessary to understand the origin 
of these differences. 

In Fig.dBwe compare the observed and the predicted fluxes 
smoothed by a Gaussian filter with a dispersion of 10 A. We nor¬ 
malized the predicted fluxes by a multiplicative factor, which 
yields the best match between observations and predictions in 
Fig . [T3] The factor was kept fixed for all wavelengths in all sub¬ 
sequent calculations. 


From the UV flux normalization we derived the stellar ra- 
dius R - 4.34 + 0.0 5 Rp, assuming a distance of 50.8 + 0.4 pc 
(Ivan Leeu wenll2007 ). This is significantly hig her than the 2.2 Rq 
reported by iPasinetti Fracassini et akl (1200 ll) or the 3.3 Rq de¬ 
rived bv iBabu ( 1977ll . On the oth er hand, ou r value agrees with 
the radius of 4.3 Rq adopted by [Rice et aP (l2004ll . Moreover, 
our derived radius yields with the rotational period and incli¬ 
nation given in Table [3] the projection of the rotational veloc¬ 
ity Vrotsinf = 53kms“', which agrees with observations (see 
TableEJ. 

The mean predicted and observed fluxes agree very well in 
the near-UV region with A > 2000 A (see Fig. (IB. However, 
the mean predicted flux is slightly lower than the observed one 
in the region 1250 - 1500 A. The opposite is true in the region 
between 1700 - 1900 A. On the other hand, the amplitude of the 
flux variations agrees well with observations in the whole UV 
domain except for the region of 2000 - 2300 A. 
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Fig. 13. Comparison of the predicted flux (mean and its variation, blue) with corresponding observed quantities (black). The pre¬ 
dicted and observed fluxes (lUE) were smoothed by a Gaussian filter with a dispersion of 10 A. Overplotted is the flux calculated 
for a slightly supersolar chemical composition with eMg = -4.1, esi = -3.7, ecr = -5.9, emii = -5.7, and epe = -4.2. 
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Fig. 14. Predicted and observed (lUE) flux in selected phases. Predicted fluxes were smoothed by a Gaussian filter with a dispersion 
of 1.3 A. Individual strong lines and iron line blends are identified. 
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Fig. 15. Comparison of the amplitudes of the observed photo¬ 
metric variability and the variability predicted by the physical 
model. The real variability is in most cases weaker with some 
exception when it is strengthened by the role of the secondary 
spot at the phase 0.11. 



1 2 3 4 5 6 7 8 


7^; [tim-'l 

Fig. 16. Dependence of the phase shift between the observed and 
predicted light curves. Note the well documented jump in the 
shift in the region of the Balmer jump. See also the shapes of 
observed and predicted light curves in Fig.fTOl 


This is further shown in Fig. [14] where we compare predicted 
and observed fluxes for selected phases. The fluxes agree very 
well in the near-UV region with A > 1900 A. The shape of the 
pseudo-continuum (formed by thousands of metallic lines) and 
the strongest line blends of mostly chromium and iron agree very 
well. The only exception is the Mg ii doublet at 2800 A, whose 
line depth is overestimated by the spectrum synthesis. 

The fluxes in the far-UV region agree less well. The pre¬ 
dicted flux is underestimated in the far UV region with A < 
1500 A with an exception of the phase around (p = 0.7. 
Moreover, the observations show a stronger line absorption in 
the region 1700 - 1800 A, which is not explained by the models. 

The fact that the observed and predicted fluxes most strongly 
disagree in the region with A < 1500 A, where the silicon 
strongly influences the flux distribution (see Fig. |7]l points to a 
possible problem with the silicon abundance map. This might 
be connected with the relatively few spectra that were used 
to derive the maps (13). This might also explain the shift be¬ 
tween predicted and observed UV light curves in the bands with 
A < 1500 A (see Fig. fTTI i. On the other hand, the stronger line 
absorption in the region 1700 - 1800 A is connected with miss¬ 
ing line opacity in the models, either caused by some element 
that was excluded from the mapping or by an incompleteness of 
the current line lists. 

A too strong predicted Mgii doublet at 2800 A can be 
caused by a vertical abundance distribution of magnesium. 
Thi s is a resonance doublet, whereas the 4481A line used 
by iKuschni^ (Il998ah for surface abundance mapping origi¬ 
nates from excited levels, where the lower-level energy of 
the corresponding transition is 8.9 eV. A similar effect, but 
weaker, is also visible in some Siii lines. Vertical abun- 
dance gradients are f requently found in other CP stars (e.g., 
Lfiftinger et alJ 200^ Rwbchikova et al.l 120081: iNesvacil et al.l 
2013t Bailey & Landstreel] 1201 3 I) . Moreover, the carbon line 


blend at 1335 A shows an equivalent width variability, which 
suggests an inhomogeneous surface distribution of this element. 

From Fig. [13] it follows that even the mean observed UV 
flux does not correspond to the solar chemical composition. 
Especially the observed far-UV flux is significantly lower than 


the solar flux. This flux decline can presumably not be mod¬ 
elled without reliably knowing the surface abundance distribu¬ 
tion. This makes the temperature determination of CP stars from 
the UV spectra problematic. 


8. Detailed comparison of observed and predicted 
light curves 

We conclude that the predicted light curves explain the nature of 
the photometric variability of 0Aur, but the correspondence is 
not excellent. The physical model, for example, does not agree 
with all results obtained by the phenomenological model de¬ 
scribed and discussed in Sect.[3] especially the incidence and 
properties of the second photometric spot centred at the phase 
0 . 11 . 

Although we see a large scatter in the ratio of the observed 
and predicted amplitudes of light variations (see Fig. [IS, the real 
light variability of the star is weaker than predicted. This might 
be because some important sources of photometric variability 
were neglected, because including additional variability sources 
results in a less important role of the remaining sources. The 
other possibility is that the variability caused by silicon is too 
large. Nevertheless, there are some interesting exceptions - the 
relatively strong enhancement of the observed amplitude in the 
photometric band u, where some mechanisms connected with 
the nearby Balmer jump might be important. 

The phase shift between the observed and predicted light 
curves in Fig. [16] provides additional information about the 
agreement between observed and predicted light curves. The 
shift in the optical region decreases with increasing frequency, 
and in the region of the Balmer jump, it suddenly switches to a 
large positive shift. The phenomenological model explains this 
phenomenon by strong changes in the amplitude of the sec¬ 
ondary spot, which reaches its highest brightness before the 
Balmer jump and almost vanishes after the jump. Trends in the 
shift in other parts of the figure can be also explained by the 
wavelength dependency of the contrast of the secondary spot. 

The predicted light curves are based on the horizon¬ 
tal elemental abundance variations derived from observations. 
However, CP stars also show vertical abundance gradients. 
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which in 0Aur are probably evidenced in some UV line pro¬ 
files. It is possible that the emerging flux at different wavelengths 
originates in different vertical regions. Consequently, the verti¬ 
cal abundance gradients in the atmosphere may in principle af¬ 
fect the light curves, and neglecting vertical abundance gradients 
may be one of the reasons for the disagreement between theoret¬ 
ical and observed light curves. 

9. Conclusions 

We successfully simulated the UV and visual flux variability of 
the chemically peculiar star 0Aur. We assumed that the vari¬ 
ability is caused by the inhomogeneous surface distribution of 
chemical elements. We used model atmospheres calculated with 
appropria t e abund ances corresponding to the Doppler maps of 
iKuschni^ (Il998al) to predict the light variability. 

We also used our phenomenological model to improve the 
rotational period of the star and to model the light curves in the 
spectral region 1300 - 5500 A. The phenomenological model 
assumes two photometric spots with distinct optical character¬ 
istics. We determined the geometry of the magnetic field and 
showed that there is no obvious connection between the field 
and photometric spots on the stellar surface. 

The numerical simulations based on the Doppler maps and 
model atmospheres explain most of the variability of 6Aur. 
Heavier elements influence the emerging flux through flux re¬ 
distribution mainly from the far-UV to near-UV and visible re¬ 
gions. The bound-free transitions of silicon and the bound-bound 
transitions of chromium, iron, and to a lesser extent manganese 
mostly cause the flux redistribution and variability in 6 Aur. As 
a result of the flux redistribution and inhomogeneous elemental 
surface distribution, the individual surface elements display dif¬ 
ferent brightnesses in individual spectral regions. The total ob¬ 
served flux, which is given by integrating over the whole stellar 
surface, is modulated by the stellar rotation. Resulting simulated 
light curves were compared with synthetic colours derived from 
lUE spectroscopy and with observed variability in bands of the 
Stromgren photometric system. 

The inhomogeneous surface distribution of silicon, 
chromium, manganese, and iron is able to explain most of the 
observed UV and visible flux variations. We reproduced the 
anti-phase behaviour of the light curves in the far-UV and visible 
regions. The light curve amplitude is about a few hundredths of 
a magnitude in the visual domain, while it reaches nearly 1 mag 
in the UV. There remains some minor disagreement between the 
predicted and observed light curves, particularly in the u colour 
of the Stromgren photometric system or in the far-UV band at 
1300 A. 

Our models reproduce the observed spectral energy distribu¬ 
tion in the near-UV region in detail. The absorption features that 
resemble strong lines in this region are in fact mostly blends of 
numerous iron and chromium lines. The spectral energy distri¬ 
bution in the far-UV region is only well reproduced around the 
phase of If - 0.7. There are minor differences in other phases, 
which may be connected with a coarser phase coverage of the 
spectra used to prepare the silicon abundance maps. Other dis¬ 
crepancies in the far-UV may be attributed to some missing line 
opacity in the models. 

The UV spectral analysis is very well suited for detecting 
vertical abundance gradients. The UV region comprises many 
strong lines that originate from states with low excitation po¬ 
tential, which provides information about the upper parts of the 
atmosphere. In combination with an optical analysis (on which 


the Doppler maps are based), which also yields information from 
the lower atmospheric parts, we were able to detect the vertical 
abundance gradient of magnesium and possibly silicon. 

Our study provides additional evidence that the light vari¬ 
ability of chemically peculiar stars is in most cases caused by 
the inhomogeneous surface distribution of individual chemical 
elements and flux redistribution and is modulated by the stellar 
rotation. For cooler chemically peculiar stars in addition to sil¬ 
icon and iron, other iron-peak elements also influence the light 
variability, in our case, mostly chromium. The comparison of the 
observed and predicted spectral energy distribution was able to 
reveal all missing opacities in the model atmospheres. 
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